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SUMMARY 

Newman Crack - closure model and the relevant crack growth program 
were applied to the analysis of crack growth under constant amplitude 
and aircraft spectrum loading on a number of aluminum alloy materials. 
The analysis was performed for available test data of 221 9-T851, 2024-T3, 
2024-T351, 7075-T651, 2324-T39 and 7150-T651 aluminum materials. 

The results showed that the constraint factor is a significant 
factor in the method. The determination of the constraint factor is 
discussed in this paper. 

For constant amplitude loading, satisfactory crack growth lives 
could be predicted. For the above aluminum specimens, the ratio of 
predicted to experimental lives, Np/Nt, ranged from 0.74 to 1. 36. The 
mean value of Np/Nt was 0.97. For a specified complex spectrum loading, 
predicted crack growth lives are not in very good agreement with the 
test data. Further effort is needed to correctly simulate the 
transition between plane strain and plane stress conditions, existing 
near the crack tip. 


INTRODUCTION 

Crack growth analysis is one of the distinct elements of damage 
tolerance design for aircraft structure. 

Prediction of crack growth under aircraft spectrum loading is very 
difficult because of retardation due to overloads, acceleration due to 
underloads which are followed by an overload and reduction of 
retardation caused by underloads following an overload. 

There are a number of crack growth models. Table 1 gives the 

comparison of the typical models (ref. 1). 

In these models, crack closure models are mathematically more 
complicated, they offer the potential to analytically predict many more 
characteristics than yield zone models and have received more and more 
attention. Newman’s crack- closure model is one of the most recent 
works on this topics. 
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This paper summarizes the analytical results of crack growth under 
constant amplitude and aircraft spectrum loading on a number of 
aluminum alloy materials. Newman crack - closure model and the relevant 
crack growth program were appl ied to the analysis ( ref. 2 ). The purpose 
of this analysis is to compare Newman’s method with test and further to 
understand its application to aircraft structural damage tolerance 

design. 


CRACK - CLOSURE MODEL 


The crack - closure model developed by Newman was based on a 
modification to the Dugdale model describing the plastic-zone ahead of 
a crack. It is assumed that the crack could be fully or partially 
closed at positive minimum stresses due to residual plas ic 

deformations left in the wake of an advancing crack. Upon reloading, 
the crack surfaces gradually separate until the crack is fully open at 
an applied stress equal to the * crack opening stress *. Figure 1 shows 

a schematic of the model. 

The crack opening stresses (Sop) were calculated from the contact 
stress at minimun load. The contact stresses were calculated from 

displacement compatibility equations with constraints added 

A constraint factor a was introduced to account for the effect of 
the state of stress on the plastic-zone size. The material is assume 
to yield in tension when the stress is a a Q and in compression when the 
stress is - or 0 .The flow stress <r 0 is taken to be the average 
between the yield stress and the ultimate tensile strength. Ideal plane 
stress or plane strain conditions are simulated with a = 1 or 3, 
respectively. In general, the constraint factor will lie between the 
ideal limits and may additionally vary with crack length. 

Newman’s crack growth rate equation is as follows: 
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The coefficients C x to C„ were determined from the best fit to 
constant amplitude experimental data. For 2219-T851 aluminum plate, 
the coefficients are ( ref. 2): 

Ci * 1. 764 X 10' 10 
C 2 = 3. 18 

C» = 2. 97 Mpa m 1/a ,-v 

C 4 = 0. 8 W 

C B = 77 Mpa m 1/a 


TEST DATA 

The analysis was performed for the available test data ( ref. 3 and 
fe .4 ). Crack growth tests were all for center cracked tension ( CCT) 
specimens. For constant amplitude loading, the test data on 2219- T851 
aluminum material are shown in Table 2 ; the test data on 2024- T3 
2024-T36! 7075-T651, 2324-T39 and 7160- T661 aluminum mateMall are 

n n n n e i e 3 * Fot spectrum loading, the test data on 2024- T3 and 
2024-T361 aluminum material are shown in Table 4. 


CRACK GROWTH RATE DATA 


To obtain the crack growth rate against AK eff relationship, 
the following methods can be used: 

1) Newman’ s crack growth rate equation ( see equation (1) ), which was 
employed fot 2219-T86I opeelmeno in this popor. 

"'“ 0d ’ wl,, ‘ "* employ ' d f " specimens except 

it 19-1861 in this paper. 

da 

dN 


The ... versus K m „ relationship was obtained from the test and 


further was transformed into ^ versus AK Bff using Newman’s crack 
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opening stress equation ( ref. 6). 


RESULTS AND DISCUSSION 

The crack growth lives from the initial crack length to the final 
crack length were calculated with various constraint factors. 

1) Results for constant amplitude loading 

The effect of the constraint factor on the ratio of predicted to 
experimental lives for various materials is summarized in Figure 2. 
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Figure 2 shows that there is a value of a, which gives the best 
comparison with calculated results for a material of a given thickness 
and stress ratio. 

For 2219-T851, Reference 2 pointed out that the equations ( 1) 
through (5) with a = 2. 3 would give a good correlation under constant 
amplitude loading. But in our analysis ( see Table 2 ), it was found 
that at the intermediate stress level ( S m »ic = 138 MPa ), the crack 
front is in a state of stress transition. As a consequence a = 1. 9 is 
found to give a better correlation. At the higher or lower stress level, 
the prediction could not best fit the test data. It seems necessary to 
fit a new set of coefficients and constraint factor. The tendency is as 
follows; at the lower stress level ( S max = 55 MPa ) - plane strain 
conditions should prevail and a should vary from 2.3 to 3. 0 as stress 
ratio R decreases; at the higher stress level ( S ma x = 276 MPa ) - plane 
stress conditions should exist and ct should approach 1.0 at positive R. 

The calculated lives are not very sensitive to a when using the 
table look - up method. For all other aluminum materials except 2219- 
T851, with R = 0. 1 ( = 124. 2, 82. 8 MPa ) , a = 1. 0 ( t < 3. 175 mm ) 

and 2. 7 ( 6. 35 mm < t < 12. 7 mm ) give the best results respectively. 

Comparison of experimental and predited crack growth lives obtained 
with the optimum a for all specimens is shown in Table 3. Figures 3 
through 5 show the crack growth analysis curves and test curves for 
three aluminum materials. 

The ratio of predicted to experimental crack growth lives, Np/Nt, 
ranged from 0.74 to 1.36 for all 29 constant amplitude data. The mean 
value of Np/Nt was 0. 97. 

2) Results for spectrum loading 

The original crack growth program only allowed input loading with 
a specified periodic load sequence. It has been modified to be able to 
analyze typical transport spectra, ‘ which have a random sequence in a 
block included thousands of flights. 

For this spectrum loading, predicted crack growth lives are not in 
very good agreement with the test data. Comparison of experimental and 
predicted crack growth lives for 2024 aluminum alloy under transport 
spectrum loading is shown in Table 4. 

The predictions made under spectrum loading were more sensitive to 
the constraint factor. Therefore, it is very important to correctly 
simulate the stress state, either plane stress or plane strain, 
existing near the crack tip. 
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CONCLUSION 


Newman crack - closure model was applied to the analysis of crack 
growth a Under simple and complex load histories for a numbed of aluminum 

The conclusions are drawn as follows; 

The key to applying the Newman crack - closure model is to 

determine a proper constraint factor a. There is a value of a 
which gives the best calculated result for a material of a given 

adpVf 88 and l ttB8a fat10, U n,eans that Effect judgement and 
adequate assumptions in the state of stress have an important effect on 
the accuracy of analysis. 

l0adi " g ’ i0 ge ” eral - satisfactory crack 
£' ’ hL Tf "’“ a once tie proper constraint factor a 

' te ” The ratio of predicted to enper taental crack growth 

Tle^enn^^l' '“J 8 ® 11 °’ 74 to L 3 ® fof 29 constant amplitude data. 

The mean value of Np/Nt was 0. 97. 

th P . P f ° r 8P : CtrtI “ laading ’ the analytical complications increase and 
the information on the transition between plane strain and plane stress 

of n ede f I 16 Constralnt factor a be a variable and a function 

f the plastic zone to thickness ratio. At short crack lengths and low 

i) 1 hop 8 t! V * B> P ] asUc zone size ia 81Ila11 compared to material thickness, 

2 and I f coitions may prevail and a may be in the range between 
„ f 3 ; ^s the crack grows longer, a will drop steadily until at some 

ger length where plastic zone size is more likely close to the 
thickness, plane stress conditions may prevail and hence a nears 1 
d«f i fjing winch could be pointed out is that the material input 
da a should cover the complete relevant range of AK and R values when 

nrn!f i! gUe " ^ Uble 88 ^ inpUt t0 tie crack growth 

program. ° 

Note that crack - cloaure model, urhich performs a cycle by cycle 

almHfl 10 ?’- J" tm - ‘ co " 8Umi ”S- 11 iB hopeful to do acme 
simplification for engineering practice. 

The following are some recommendations for the program; 

As the calculated crack opening stress was held constant in the 

program while the crack growth over the length AC*, the crack growth 

accumulation procedure may be simplified. 

’ e^ivalent * crack opening stress might be used as early as 
possible in order to reduce the computer time. 
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Table 1. Comparison of Typical Crack Growth Models 


TYPE 

FEATURE PARAMETERS 

CONTRIBUTORS 


[1] Plastic Yield Zone 
Mode 1 

Effective Stress 
Intensity Range 

Wheeler (1970) 

Wi 1 lenborg (1971) 
Gallagher (1974) 
Chang (1981) 
Johnson (1981) 


[2] Crack Closure 
Mode 1 

Crack Opening Stress 

Elber (1969) 

Newman (F in i te element 
(Contact stress 
Bell (AFFDL) (1974) 
Maarse (1977) 

Matsuoka (1976) 

De Koning (1980) 

1974) 

1981) 

[ 3] Modi f ied Model 
Based on One or 
Both of above Types 

One or Both of [1] 
and [2] 

Huang (1980) 
He (1980) 
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Table 2. Comparison of Experimental and Predicted Crack Growth 

Lives for 2219-T851 Aluminum Alloy 


Smax 

MPa 

R 

at 

mm 

flf 

mm 

a =1. 0 

Np/Nt 
a =1. 9 

a =2. 3 

a =2. 7 

55 

-0. 1 

3. 91 

18. 7 



1. 62 

1. 36 

55 

0. 0 

3. 55 

52. 4 



0. 99 

0. 83 

55 

0. 3 

3. 87 

6. 9 



0. 76 

0. 64 

138 

-0. 3 

4. 13 

45. 7 


1. 11 



138 

-0. 1 

3. 85 

43. 0 


0. 99 



138 

0. 01 

3. 87 

50. 0 


0. 99 

0. 83 


138 

0. 01 

4. 00 

42. 0 


1. 04 



138 

0. 01 

6. 67 

38. 3 


0. 95 



138 

0. 2 

4. 57 

48. 0 


1. 22 



138 

0. 3 

4. 06 

43. 6 


0. 80 



138 

0. 7 

3. 94 

12. 1 

1. 12 

0. 87 



276 

-0. 1 

3. 94 

21.4 



1. 15 

1. 09 

276 

0. 01 

3. 83 

14. 9 

0. 75 

0. 75 

0. 72 


276 

0. 3 

3. 94 

12. 9 

0. 74 

0. 60 

0. 58 


276 

0. 7 

3. 94 

23. 5 

0. 86 




Note: R 

- Stress 

Ratio ( Smin/Smax) 





a t 

- Half Length of 

initial crack 





a f 

- Half Length of 

final crack 





Np 

- Number 

of cycles predicted from analysis 



Nt 

- Number 

of cycle 

s from test 
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Table 3. Comparison of Experimental and Predicted Crack Growth 
Lives for Various Aluminum Mater ials under Constant Amplitude Loading 


Material 

Temper 

Smax 

(MPa) 

R 

t 

(mm) 

at 

(mm) 

Ef 

(mm) 

Np/Nt 

a 

2024 

T3 

124. 2 

0. 1 

3. 175 

7. 214 

57. 302 

0. 92 

1. 0 


T3 

124. 2 

0. 1 

1. 600 

8. 687 

59. 715 

0. 96 

1. 0 


T351 

124. 2 

0. 1 

3. 175 

10. 008 

57. 302 

0. 97 

1. 0 


T351 

124. 2 

0. 1 

6. 350 

8. 992 

56. 667 

0. 93 

2. 7 


T351 

124. 2 

0. 1 

12. 70 

8. 484 

58. 522 

0. 98 

2. 7 

7075 

T651 

82. 8 

0. 1 

3. 175 

7. 366 

60. 452 

0. 97 

1. 0 


T651 

82. 8 

0. 1 

6. 350 

6. 655 

58. 293 

0. 99 

2. 7 


T651 

82. 8 

0. 1 

12. 70 

6. 985 

58. 141 

0. 99 

2. 7 

2324 

T39 

82. 8 

0. 1 

3. 175 

14. 808 

102. 438 

0. 93 

1. 0 


T39 

82. 8 

0. 1 

7. 087 

16. 866 

114. 071 

0. 95 

2. 7 


T39 

82. 8 

0. 1 

12. 649 

25. 375 

112. 293 

0. 99 

2. 7 

7150 

T39 

82. 8 

0. 1 

3. 302 

15. 773 

111. 836 

0. 95 

1. 0 


T39 

82. 8 

0. 1 

6. 350 

15. 926 

87. 401 

0. 96 

2. 7 


T39 

82. 8 

0. 1 

12. 70 

17. 475 

117. 450 

0. 99 

2. 7 


Table 4. Comparison of Experimental i 
Lives for 2024 Aluminum Material under 

md Predicted Crack Growth 
Transport Spectrum Loading 

Temper 

t a i a f 

(mm) (mm) (mm) 

Np/Nt 

T3 

3.175 6.528 25.349 

1. 52 

T351 

6.350 5.842 60.198 

0. 30 
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2219r-I85 t . 


The number under X axis 
is the value of <x 



Figure 2. The effect of the constraint factor a on predicted lives. 
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Nornallzed crack length 



Figure 3. Crack growth comparison between prediction and test for 2219 - T851. 
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